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Radio structures of the nuclei of nearby Seyfert galaxies 
and the nature of the missing diffuse emission 
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ABSTRACT 

We present archival high spatial resolution VLA and VLBA data of the nuclei of seven 
of the nearest and brightest Seyfert galaxies in the Southern Hemisphere. At VLA 
resolution (~ 0.1 arcsec), the nucleus of the Seyfert galaxies is unresolved, with the 
exception of MCG-5-23-16 and NGC 7469 showing a core-jet structure. Three Seyfert 
nuclei are surrounded by diffuse radio emission related to star-forming regions. VLBA 
observations with parsec-scale resolution pointed out that in MRK 1239 the nucleus 
is clearly resolved in two components separated by ~ 30 pc, while the nucleus of 
NGC 3783 is unresolved. Further comparison between VLA and VLBA data of these 
two sources shows that the flux density at parsec scales is only 20% of that measured by 
the VLA. This suggests that the radio emission is not concentrated in a single central 
component, as in elliptical radio galaxies, and an additional low-surface brightness 
component must be present. A comparison between Seyfert nuclei with different radio 
spectra points out that the "presence" of undetected flux on milli-arcsecond scale is 
common in steep-spectrum objects, while in flat-spectrum objects essentially all the 
radio emission is recovered. In the steep-spectrum objects, the nature of this "missing" 
flux is likely due to non-thermal AGN-related radiation, perhaps from a jet that gets 
disrupted in Seyfert galaxies because of the denser environment of their spiral hosts. 
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1 INTRODUCTION 

Only a small fraction (~ 10%) of the population of ac- 
tive galactic nuclei (AGN) possess a powerful radio emission 
(Liaghz > 10 23 W/Hz), as found in radio galaxies/quasars 
and blazars. Seyfert galaxies are part of the "radio-quiet" 
AGN population, with radio luminosity Li.4GHz ^ iq 20 ~ 23 
W/Hz. Despite their weak radio emission, Seyfert nuclei are 
very nearby, allowing the study in detail of the radio prop- 
erties of their central engine. 

Radio observations with arcs econd resolution o f sev- 
eral Seyfert samples (see e.g. Ulvestad fc Wilson] Il984l ; 



structures reaching hundreds of kpc or even Mpc scales, 
while in Seyferts the radio emission is confined to a few 
kpc or even sub-kpc scales. 

When observed with milli-arcsecond resolution, the pc-scale 
structure of Seyfert nuclei is usually resolved i n several 



components (e.g. NGC 3079, iTrotter et ail Il99&t) 



iMorganti et al.ll 19991 : iThean et al Il2000h showed that a large 
fraction of Seyferts have resolved structures, with hint of 
jets and/or extended emission, the latter usually related 
to star-formin g regions. Se veral objects, such a s NGC 1052 
' Wrobe]||l984h. NGC 1 068 (|Ulvestad et alJll987t>. NGC 7674 



Momiian et al] l2003h . and MRK 3 l|Kukula et all Il999h 



have been found to display a radio morphology with core, 
collimated jets and hot spots, similar to those found in radio- 
loud galaxies. However, powerful radio sources have linear 
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bling a jet struct ure (e.g. NGC 4151. lUlvestad et al.lll998l : 
iNagar et al.ll200lh and sometim es with the presence of ex- 
tended emission (e.g. NGC 5793. [Hagiwara etalll2000l) . The 
comparison between arcsec and milli-arcsec radio properties 
pinpointed a frequent misalignment between pc and kpc- 
scale jets, suggesting either a change in jet ejection axis, or 
a bending due to pressu re gradients in the ambient medium 
jMiddelberg et alj|2004h . 

An intriguing characteristic shown by a large number of 
Seyfert nuclei is that the radio emission arising from their 
pc-scale structure is often much fainter than that derived 
from observations with lower resolution, even in the case the 
nucleus is unresolved. This result suggests that in Seyfert 
nuclei the radio emission is not concentrated in the cen- 
tral region, as found in powerful radio galaxie s, but it ex- 
tend s on scales of tens or hundreds pc (see e.g. ISadler et al.l 
Il995lh However, not all the Seyfert nuclei have missing flux 
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on parsec scales, as in the case of MRK 530 l|Lal et al.ll20ol) . 
indicating that the radio emission mainly arises from the 
central compact component, without evidence of extended, 
low-surface brightness features. 

In this paper, we present the results of multi-frequency 
archival VLA and/or VLB A data of a sample of some 
of the nearest and brightest Seyfert galaxies, taken from 
the infrared high spati al resolution studies co nducted by 
iPrieto et all <|20Qgf > and iReunanen et all i|2009D . For these 
sources, no complete and/or unambiguous information on 
their radio properties could be found in the literature. The 
comparison between these data with those at pc-scale res- 
olution available from the literature allows a better deter- 
mination of the physical conditions of the radio emission at 
kpc and pc scales. 

Throughout this paper, we assume the following cosmol- 
ogy: H = 71kms~ 1 Mpc" 1 , fl u = 0.27, and fi A = 0.73, in 
a flat Universe. The spectral index is defined as S(v) oc v~ a . 



2 THE SAMPLE 

The objects studie d in this pap e r are part of the sample of 
Seyfer t nuclei from lPrieto et ail (|2009l ) and IReunanen et al.l 
(2009). They consist in the nearest and brightest Seyfert 
nuclei accessible from the Southern Hemisphere. The prox- 
imity of these objects allows us to study them at scales of a 
few tens of parsec in the IR with VLT - adaptive optics and 
interferometry which deliver spatial resolutions comparable 
to that achieved in the optical with HST, and in radio with 
VLA and VLBI. In Fig. [Q we present VLT-VISIR images 
at 11.8 /im of the Seyfert studied in this paper. Details on 
VISIR observat i ons an d data reduction can be found in 
IReunanen et all (HH). The dominant IR emission comes 
from the central nucleus that is unresolved with FWHM < 
0.35 arcsec, which corresponds to linear sizes between 30 
and 190 pc, depending on the redshift. In NGC 7469 and 
NGC 7582, an additional emission from the star forming 
regions is clearly detected. 

High spatial resolution spectral energy distribution (SED) 
have been produced for these objects. For the majority 
of the Seyfert considered, multi-frequency, high-resolution 
information in the radio band has not been exploited in full. 
For this reason, we collected and analysed archival radio 
data with (sub-)arcsecond resolution in order to implement 
the information already available and to constrain the radio 
emission from these nuclei. 

The list of the Seyfert nuclei studied in this paper is in 
Table Q] We considered only datasets at frequencies and 
resolution not available in the literature, or for which no 
unambiguous information could be found. 



2.1 Radio properties of the selected sources 

2.1.1 NGC1097 

NGC 10 97 is a barred spiral galaxy hosting a Seyfert 1 
nucleus (|Storchi-Bergmann et al.1ll993t ). The radio emission 
is dominated by a ring-like structure related to star-forming 
regions also detected in optical. At the centre of the ring 



there is an unresolved (beam ~ 0.3 arcsec at 8.4 GHz) com- 
ponent jThean et al. 2000) with an inverted non-thermal 
spectrum dMorganti et all 1 19991 ). The radio core was ob- 
served by ISadler et alj (|l995l ) with the Parkes-Tidbinbilla 
interferometer (PTI) that provides a resolution <0.1 arcsec. 
No radio emission was detected, implying that the core flux 
density is <5 mjy at 2.3 GHz. 



2.1.2 MCG-5-23-16 

MCG-5-23-16 is a SO galaxy host ing a Seyfert 2 nucleus. Its 
radio properties were studied bv lUlvestad fc Wilson! (Il984l ) 
with VLA observations at 1.4 and 4.8 GHz. The radio 
structure of the nucleus is found slightly resolved (beam ~ 
0.4 arcsec at 4.8 GHz) with a spectral index a ~ 0.5. No 
evidence of extra- nuclear emission was detected. 
High spat ial resolution (bea m < 0.1 arcsec) observations 
with PTI (Sa dler et ai1 ll995) could set only an upper limit 
to the flux density of the core 5*8.4GHz <7 mjy. 



2.1.3 MRK 1239 

MRK 1239 is an early-type E/S0 galaxy hosting a Seyfert 
1.5 nucleus. It s cent ral ra dio structur e was studied by 
lUlvestad et~al1 (|l995l ) and iRush et all (|l996l ) with the 
VLA. The source appears unresolved at all frequencies 
(beam ~ 0.2 arcsec at 8.4 GHz) witho ut evidence of 
extend ed emission surrounding the nucleus. lUlvestad et alj 
|l995l ) noted a strong steepness in the radio spectrum 
between 4.8 and 8.4 GHz (a ~ 1.6), that they interpreted 
as due to flux variability (observations were not carried out 
simultaneously) . 



2.14 NGC 3783 

NGC 3783 is a barred spiral galaxy with a highly variable 
Sey fert 1 nucleus. Mult i- freque nc y radio VLA observ ations 
by lUlvestad fc Wilson! (| 19841 ) . lunger et al] l|l986l ) and 
ISchmitt et al. ~ if200l[ ) pointed out the presence of an unre- 
solved component (beam ~0.25 arcsec at 8.4 GHz) with a 
spectral index a ~ 0.5. No extended extra-nuclear emission 
was found. High spat ial resolution (bea m < 0.1 arcsec) 
observations with PTI l|Sadler et al. I l 19951 ) could set only an 
upper limit to the flux density of the core <Si.6GHz <6 mjy. 



2.1.5 NGC 5506 

NGC 5506 is a spiral galaxy with a Sey fert 2 nucleus. 
The c entra l radio structure wa s studied by lUlvestad et al.l 
j|l98ll ) and lUnger et all l| 19861 ) by VLA and MERLIN ob- 
servations. It shows a diffuse bubble-like feature extending 
to the north-west of t he unresolved nucleus. At 5 GHz 
IWehrle fc Morris! l| 19871) found a low-surface brightness halo 
enshrouding the central structures. 

Multi- frequency observa tions of the nucleus with the 
PTI (|Sadler et all Il995h showed that it has a convex 
spectrum. Pc-scale VLBI observations at 18, 6 and 3.6 cm 
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Figure 1. VLT-VISIR images at 11.8 /im of the Seyfert nuclei studied in this paper. Images adapted from Reuna nen et al . (2009). 
Contours are at 3cr, 5, 7, 11 and 19 levels. The triangle in NGC 7582 represents the ionization cone. 



jMiddelberg et all [2004) could resolve the nucleus in three 
main components roughly aligned in the east-west direction 
and extending within an area of about 50 mas (~ 6 pc). 
The flattish spectrum displayed by the brightest component 
suggests that it is the true source core. 



2.1.6 NGC 7469 

NGC 7469 is a spiral galaxy with a Seyfert 1 nucleus 
embedded in a ring of starburst activity. At radio frequen- 
cies, its structure shows an unresolved centr al component 
surro unded by a ring of star-forming regions (|Wilson et al.l 
1991). MERLIN observations with high spatial resolution 
(beam ~ 0.05 arcsec) showed that the nucleus is clearly 
resolved in two components with a core-jet structure elon- 
gated in the east-west direction. VLBI observations (beam 
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Figure 2. VLA images at 8.4 GHz {top) and at 14.9 GHz {bottom) 
of the central region of NGC 1097. On each image we provide the 
observing frequency; the restoring beam, plotted on the bottom 
left corner; the peak flux density in mJy/beam; the first contour 
intensity {f.c. in mJy/beam), that is 3 times the off-source noise 
level; contour levels increase of a factor 2. 
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Figure 3. VLA image at 8.4 GHz of the central region of MCG- 
5-23-16. On the image we provide the observing frequency; the 
restoring beam, plotted on the bottom left corner; the peak 
flux density in mJy/beam; the first contour intensity {f.c. in 
mJy/beam), that is 3 times the off-source noise level; contour 
levels increase of a factor 2. 



Table 1. The Seyfert sample. Column 1: source name; col. 2: 
type; col. 3: redshift; col. 4: luminosity distance; col. 5: scale. 



Source 


Type 


z 


D L 


kpc/" 


NGC 1097 


1 


0.00424 


18.0 


0.086 


MCG-5-23-16 


2 


0.008486 


36.1 


0.172 


MRK 1239 


1.5 


0.02883 


124.5 


0.570 


NGC 3783 


1 


0.009730 


41.4 


0.197 


NGC 5506 


2 


0.006181 


26.2 


0.126 


NGC 7469 


1 


0.016317 


69.8 


0.328 


NGC 7582 


2 


0.005254 


22.3 


0.107 



~ 10 mas ) at 1.6 GHz (|Lonsdale et al.ll2003h could resolve 
the core-jet structure in five different components, lying in 
an east-west line as found in the lower resolution MERLIN 
image, and contained in an area of about 168 mas (~ 55 pc). 



spati al resolution (beam < 0.1 arcsec) observations with 
PTI (|Sadler et all Il995h could set only an upper limit to 
the flux density of the core Ssaghz < 5 mjy. 



2.1.7 NGC 7582 

NGC 7582 is a barred spiral galaxy hosting a Seyfert 2 
nucleus surrounded by a star forming ring. Its radio struc- 
ture shows an unresolved central component embedded in 
a diffuse emission related to the star-forming regio n also 
detected in the optical (|Ulvestad fc Wilsoni Tl984h . The 
extended emission is elongate d in the NW-SE dir ection, 
with a spectral index a ~ 0.7 (|Morganti et al.lll999l V High 



3 NEW RADIO DATA 

Archival VLA/ VLB A radio data of the aforementioned 
sample were analysed in order to implement the information 
available in the literature. To achieve an adequate spatial 
resolution to resolve the nuclear structure, we considered 
observations carried out when the VLA was in one of 
the extended configuration (array A or B) at frequencies 
ranging between 1.4 and 14.9 GHz, when available. We 
reduced datasets either at frequencies lacking information 
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Figure 4. VLA images at 8.4 GHz (top), and 14.9 GHz (cen- 
tre), and VLB A image at 1.6 GHz (bottom) of the central re- 
gion of MRK 1239. On each image we provide the observing fre- 
quency; the restoring beam, plotted on the bottom left corner; the 
peak flux density in mJy/beam; the first contour intensity (f.c. 
in mJy/beam), that is 3 times the off-source noise level; contour 
levels increase of a factor 2. 



Figure 5. VLA images at 4.8 GHz (top), 8.4 GHz (centre), and 
14.9 GHz (bottom) of the central region of NGC 5506. On each im- 
age we provide the observing frequency; the restoring beam, plot- 
ted on the bottom left corner; the peak flux density in mJy/beam; 
the first contour intensity (f.c. in mJy/beam), that is 3 times the 
off-source noise level; contour levels increase of a factor 2. 
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or for which data presented in the literature were not found 
satisfactory. 

The data reduction was carried out following the standard 
procedures for the VLA implemented in the NRAO AIPS 
package. Images were produced after a few phase-only cali- 
bration iterations. In order to obtain accurate flux density 
at 1.4 GHz, it was necessary to image several confusing 
sources falling within the primary beam. Uncertainties 
in the determination of flux density are dominated by 
amplitude calibration errors, which are between 3% and 
5%, being worse at 14.9 GHz. The rms noise level (la) on 
the image plane is usually below 0.1 mjy, being irrelevant 
for our targets, with the exceptions of NGC 1097 and 
MCG-5-23-16 where it is comparable to the amplitude 
calibration errors. 

In the case of NGC 1097, MRK 1239 and NGC 3783, for 
which no information on their pc-scale structure could 
be found in the literature, we reduced archival VLBA 
data. No archival VLBA data were found for the sources 
MCG-5-23-16 and NGC 7582. The Seyfert nuclei considered 
in this paper are not bright enough to give detectable 
fringes. VLBI fringe fitting was done on a bright source 
lying close to the target, and off-beam phase referencing 
was performed to calibrate the phases of the target source. 
In the case of NGC 1097, no signal of the source was 
detected. For MRK 1239 and NGC 3783, the structure was 
visible from the first image with a 4a noise level. Then 
phase self-calibration with a solution interval of 30s was 
performed using the CLEAN component model. After a 
few iterations we use natural weights to pick up possible 
extended emission. No amplitude calibration was performed 
due to the weakness of the targets. 

4 RESULTS 

4.1 Radio images 

Full resolution VLA images of the sources discussed in this 
paper are presented in Figs. [2] to [7] In the case of NGC 7469 
and NGC 7582, where an extended diffuse emission is 
present around the central nucleus, an image without the 
shortest baselines was produced in order to better describe 
the compact components without the contamination from 
the low-surface brightness extended emission (Figs. [Tip and 
Eb). For MRK 1239 and NGC 3783, VLBA images at 1.6 
GHz are presented in Figs.[4j: and [8] respectively. 
Flux density and angular size were measured by means 
of the task JMFIT, which performs a Gaussian fit on 
the image plane. In case of extended structures, the flux 
density was derived by TVSTAT, while the source size was 
measured from the lowest contour on the image plane. The 
nucleus of these Seyfert galaxies is unresolved or marginally 
resolved in VLA observations. We consider marginally 
resolved those nuclei whose largest angular size (LAS) is 
between 0.5 and 1 beam size at the best resolution, and we 
term unresolved those whose LAS is smaller than half of 
the beam. 

Source parameters, together with information on the 
observations, are reported in Tabled In Table 3, we present 
the observational parameters of the source components. 



Table 3. Radio properties of the source components. Col. 1: 
source name; col. 2: source component, the name refers to the 
label as reported in Figs. from[2]to[8] col. 3: observing frequency; 
col. 4: flux density; cols. 5 and 6: deconvolved major and minor 
axis, a: values derived on VLBA images. 



Source 


Comp. 


Freq. 


Flux 


^max 








GHz 


mjy 


arcsec 


arcsec 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


TVPr 1 1 HQ7 
in vj v_ ivjy ( 


n 

V j 


R A 

O.T 


1 n 

l.U 








Q 


14.9 


5 6 


<0.5 






R 1 


8 4 


1 8 


9 


0.7 




K2 


8 4 


8 


1 7 


1 3 




R3 


8 4 


6 


1 6 


8 


1V1 \y \J~ iJ-ZO - J- U 




8 A 


2 1 


<^n i 






Ij 


8 4 


5 


■\u. ± 






Q 


1 4 


60.0 


< 0.3 






Q 


5 


26.8 


<--(") i 

*^u. ± 






Q 


8 4 


15.5 


\U. 1 








14.9 


9.3 


<0.2 






A° 


1 6 


6 


0.013 


0.006 




13" 


1.6 


4.2 


0.013 


0.004 


Nnr* 378'} 


C° 


1 6 


4 6 


0.022 


0.007 


NCO 5506 




1.4 


304.0 


<1.0 






Q 


4 8 


136 








p 
V i 


R 4 

Ci.-l 


ra n 








C 


14.9 


58.0 


<0.2 






Wing 


4.8 


27 


2.1 


2.0 




Wing 


8.4 


16 


2.2 


1.7 




Halo 


4.8 


18 


4.2 


.3.0 




Halo 


8.4 


11 


3.8 


.3.2 


NGC 7469 


C 


8.4 


14.8 


<0.14 






C 


14.9 


10.6 


0.1 


<0.05 




Rl 


8.4 


0.3 


0.8 


0.4 




R2 


8.4 


0.6 


0.9 


0.7 




R3 


8.4 


0.7 


1.1 


0.8 


NGC 7582 


C 


1.9 


9.5 


0.8 


<0.1 




C 


8.4 


6.9 


0.1 


0.2 




Rl 


4.9 


3.6 


0.4 


0.3 




Rl 


8.4 


2.1 


0.3 


0.3 




R2 


4.9 


0.8 


0.3 


<0.1 




R3 


4.9 


1.1 


0.5 


0.3 



4.2 Notes on individual sources 

In this section we describe the characteristics of the Seyfert 
nuclei arising from the new images presented in this paper. 



4.2.1 NGC1097 

The new images at 8.4 and 14.9 GHz are presented in Fig. 
H and have a resolution of 0". 66x0". 25 and 1". 15x0". 45 
respectively. Fig. [2^ shows for the first time the ring 
structure at 8.4 GHz. Four star-forming regions, labelled 
Rl, R2, R3 a nd R4 in Fig. [2^,, are also present at 1.4 and 5 
GHz (see e.g. iHummel et al] [l987l. The central component 
appears unresolved in both 8.4 and 14.9 GHz images, giving 
an upper limit to its linear size of < 0".5 (<43 pc). The 
spectral index is found inverted between 1.4 and 14.9 GHz. 
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Table 2. Radio observations and source parameters for the Seyfert galaxies. Column 1: source name; col. 2: observing frequency; col. 
3: array configuration; col. 4: beam size; col. 5: observing date; col. 6: flux density; col. 7: luminosity; col. 8: largest angular size; col. 9: 
largest linear size; col. 10: morphology of the nuclear component (Un=unresolved, MR= marginally resolved; Double=two well-separated 
components); col. 11: morphology of the extended structure, when present. 



Source Freq. conf beam Obs. date Flux density Log P LAS LLS Morph. Ext 

GHz arcsec mjy W/Hz arcsec pc 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 



NGC1097 


8.4 


A 


0.66x0.25 


Oct 31 1992 


NGC 1097 


14.9 


B 


1.15x0.45 


Feb 24 2001 


MCG-5-23-16 


8.4 


A 


0.43x0.23 


Aug 21 1999 


MRK 1239 


1.4 


A 


1.54x1.25 


Jan 16 1993 


MRK 1239 


4.8 


A 


0.46x0.40 


Jul 31 1987 


MRK 1239 


8.4 


A 


0.33x0.26 


Jan 16 1993 


MRK 1239 


8.4 


A 


0.26x0.22 


Jul 15 1995 


MRK 1239 


14.9 


B 


0.97x0.42 


Dec 15 2003 


MRK 1239 


1.6 


VLBA 


0.011x0.004 


May 09 2005 


NGC 3783 


1.6 


VLBA 


0.028x0.008 


Nov 30 1999 


NGC 5506 


1.4 


A 


1.71x1.24 


May 06 2002 


NGC 5506 


4.8 


A 


0.48x0.37 


May 06 2002 


NGC 5506 


8.4 


A 


0.27x0.20 


Nov 14 2000 


NGC 5506 


14.9 


B 


0.72x0.43 


Dec 05 1991 


NGC 7469 


8.4 


A 


0.22x0.19 


Feb 05 2006 


NGC 7469 


14.9 


A 


0.14x0.10 


Sep 08 1999 


NGC 7582 


4.8 


A 


1.25x0.32 


Dec 12 2000 


NGC 7582 


8.4 


A 


0.75x0.19 


Dec 12 2000 



30.0 


21.04 


18 


1550 


I .hi 


Ring 


5.0 


20.30 


<0.5 


< 48 


Un 


No 


2.6 


20.60 


1.2 


206 


MR 


No 


60.0 


23.04 


<0.3 


<171 


Un 


No 


27.0 


22.60 


<0.1 


<57 


I in 


No 


15.0 


22.44 


<0.1 


<57 


I in 


No 


14.5 


22.44 


<0.1 


<57 


Un 


No 


9.0 


22.22 


<0.2 


114 


Un 


No 


10.2 


22.26 


0.070 


40 


Double 


No 


4.6 
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4.2.2 MCG-5-23-16 

The VLA image at 8.4 GHz (Fig. with a resolution of 
0". 43x0". 23, is presented for the first time. The source is 
resolved in the north direction, suggesting the presence of a 
jet-like structure. The total radio spectrum shows a slight 
steepening at 8.4 GHz with a spectral index a = 0.8 ± 0.1. 



4.2.3 MRK1239 

New 14.9-GHz VLA and 1.6-GHz VLBA images are 
presented in Figs. |4p,c with resolution of 0". 97x0". 42 
and 0". 011x0". 004 respectively. The central structure is 
unresolved in our VLA images at all frequencies, giving 
an upper limit to the linear size of < 57 pc (Figs. |4^,,b). 
The spectral index we obtain is a ~ 0.9 ± 0.1 without 
any ev idence of the steepening reported bv lUlvestad et all 
<|l995l) (see Section 2.1.3), an d in a better agreement 
with what found by iRush et al. I (Il996t). It is worth noting 



that data reported by Ru sh et al. I ll 19961 ) refer to low 



resolution observations carried out when the VLA was in 
D-configuration. The good agreement between these values 
and our measurements suggests that no extended emission 
on kpc scales is present. When observed with pc-scale 
resolution, the nucleus of MRK 1239 at 1.6 GHz is resolved 
in two components separated by ~ 50 mas ( ~ 30 pc), with 
position angle of 40° (Fig. The lack of multi-frequency 
observations with similar resolution does not allow us to 
study the spectral index of these components. 
The flux density measured on the VLBA image at 1.6 GHz 
is 10.2 mjy, i.e. only 20% of the VLA flux density at 1.6 
GHz, obtained re-scaling the 1.4 GHz VLA flux density 



with the spectral index computed between 1.4 and 4.8 
GHz (a = 0.6). This indicates that almost 80% of the flux 
density measured with the VLA is missing in the VLBA 
image. 



4.2.4 NGC3783 

The VLBA image at 1.6 GHz, with a resolution of 
0". 028x0" .008 (Fig. O, is presented for the first time. In 
this image, the nucleus is unresolved with a linear size < 4 
pc. The total flux density measured on the VLBA image is 
4.6 mjy, that represents only a 20% of the flux density of the 
unres olved component in VLA observations (lljnger et al.l 
Il98rj ). This indicates that the majority (~80%) of the flux 
measured with the VLA is missing in the VLBA image. 



4.2.5 NGC5506 

New VLA images at 4.8 and 8.4 GHz with resolution 
of 0".48x0".37 and 0". 27x0". 20 respectively, are shown 
in Fig. [5^,b. The source displays an unresolved central 
component with a diffuse wing-like emission extending 
mainly to the north-west and east part of the nucleus. 
The high dynamic range, namely the ratio between the 
peak flux density and lcr noise level, of our images allows 
us to detect also at 8.4 GHz the extended low-surface 
brightness halo (Figs. [5^,b). This has a diameter of ~ 
2". 75 (~ 350 pc) enshrouding the central features. The 
nucleus is unresolved in our VLA images, giving an upper 
limit of <0".08 (<10 pc). The nucleus accounts for the 
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majority (75%) of the radio emission. Its spectral index 
is a ~ 0.8 ± 0.1. The low-surface brightness halo and 
the extended wing-like structures have a slightly steeper 
spectral index a — 0.9 ± 0.1. 

A comparison between simultaneous MERLIN and EVN 
observations at 18 and 6 cm clearly indicates that almost 
43% of the flux density detected by MERLIN c annot be 
recovered in EVN images (|Middelberg et ahl 120041 ). In fact, 
a diffuse emission, detected by MERLIN observations, is 
not seen by the EVN, as the EVN is insensitive to struc- 
tures larger than 35 and 11 mas at 18 and 6 cm, respectively. 



4.2.6 NGC 7469 

We produced new VLA image at 8.4 GHz with a resolution 
of 0". 2x0". 18 (Fig. |6p) obtained without the shortest 
(< 35 kA) baselines in order to pinpoint the compact 
central structure and to reduce the contamination from 
the surrounding ring. The radio emission is dominated by 
an unresolved central component surrounded by a ring 
of star-forming regions with a diameter of about 3 .7 (~ 
1.2 kpc). In our new image at 14.9 GHz (Fig. [6};), the 
nucleus is elongated in the east- west direction. The spectral 
index is a ~ 0.5 ±0.1. The el ongation is in agree ment 
with t hat shown by MERL IN (jAlberdi et al.l l200rj) and 
VLBI (|Lonsdale et al.l I2003T I images where the nucleus 
clearly shows a core-jet structure elongated in the same 
direction. In the diffuse extranuclear emission, three clumps 
of star-forming regions, labelled Rl, R2, and R3 in Figs. 
[6^,b, are clearly visible at 8.4 GHz. 

A comparison between MERLIN and VLBA images clearly 
indicates that almost 50% o f the flux density detected by 
MERLIN (FWHM ~ 0".15. iThean et al.ll200ll) is missing 
in th e VLBA image at the same frequency (Lonsda le et al.l 
l2003h . 
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Figure 6. VLA images of the central region of NGC 7469 at 
8.4 GHz (top), without the shortest baselines (centre), and at 
14.9 GHz (bottom). On each image we provide the observing fre- 
quency; the restoring beam, plotted on the bottom left corner; the 
peak flux density in mJy/beam; the first contour intensity (/. c. 
in mJy/beam), that is 3 times the off-source noise level; contour 
levels increase of a factor 2. 



4.2.7 NGC7582 

We produced new VLA image at 4.8 GHz with a resolution 
of 0". 94x0". 24 (Fig. [7p) obtained without the shortest 
(< 45 kA) baselines in order to pinpoint the compact 
central structure and to reduce the contamination from 
the surrounding extended emission. The radio emission is 
dominated by an unresolved central component surrounded 
by diffuse emission related to star- forming regions (Fig. [7}, 
also detected in the optical. The nucleus, with a spectral 
index a~0.6±0.1, is unresolved at all frequencies, giving 
an upper limit of about 0".4 (<40 pc) to its linear size. The 
extended emission has a maximum diameter of about 9 (~ 
950 pc) and a spectral index a ~ 1 .0 ± 0.1, that is steepe r 
than the value (a ~ 0.7) found by iMorganti et all (|1999T) . 
The discrepancy may arise from the different resolution of 
the observations taken into account. Three compact regions, 
labelled Rl, R2, R3 in Figs. [7£i and[7p, are present in the 
diffuse emission at 4.8 GHz. For these clumps, near-infrared 
counterpart was found by Fernandez-Ontiveros (2009, in 
preparation) suggesting that they are star-forming regions. 
Component Rl is visible also at 8.4 GHz, with a flux 
density of ~2.1 mjy, and a spectral index a ~ 0.9 ± 0.2, 
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Figure 7. VLA images of the central region of NGC 7582 at 
4.8 GHz (top), without the shortest baselines (centre), and at 
8.4 GHz (bottom). On each image we provide the observing fre- 
quency; the restoring beam, plotted on the bottom left corner; the 
peak flux density in mJy/beam; the first contour intensity (/. c. 
in mJy/beam), that is 3 times the off-source noise level; contour 
levels increase of a factor 2. 
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Figure 8. VLBA image at 1.6 GHz of the central region of 
NGC 3783. On the image we provide the observing frequency; 
the restoring beam, plotted on the bottom left corner; the peak 
flux density in mJy/beam; the first contour intensity (/. c. in 
mJy/beam), that is 3 times the off-source noise level; contour 
levels increase of a factor 2. 



that is a little flatter than the spectral index of the whole 
extended emission. 



5 DISCUSSION 

5.1 The nuclear radio emission 

In general, the nucleus of these Seyfert galaxies is unresolved 
at the resolution provided by the VLA in its larger configu- 
rations (A and B arrays), giving an upper limit to the linear 
size between about 20 and 120 pc. Only MCG-5-23-16 and 
NGC 7469 are marginally resolved at the highest resolution 
(FWHM ~ 0.2 arcsec). In both cases the central component 
is elongated suggesting a core-jet morphology. 
For all the sources we find that the spectral index of the 
nucleus ranges between 0.5 and 0.9, with the only exception 
of NGC 1097 where the spectrum is inverted. 
When observed with pc-scale resolution, most sources dis- 
play a resolved st ructure. This is the case of MRK 1239 (Fig. 
lib), NGC 5506 dMiddelberg et al.l [20o3 ), and NGC 7469 
( Lonsdale et~aHl2003| >. MRK 1239 shows a well defined dou- 
ble morphology with a total linear size of 40 pc. The lack of 
mult i- frequency information does not allow us to study the 
spectral index distribution. NGC 5506 and NGC 7469 are 
characterised by several compact components resembling a 
core-jet structure. In NGC 5506, the availability of multi- 
frequency observatio ns singles out one of th e components 
with a fiat spectrum (jMiddelberg et alj|2004h . 
NGC 3783 is still unresolved at pc-scale resolution (Fig. [BJ, 
and an upper limit to the linear size < 4 pc is estimated. 
For the sources with VLBA information, MRK 1239, 
NGC 3783, NGC 5506, and NGC 7469, an estimate of the 
brightness temperature Tb, the equipartition magnetic field 
H ecl and minimum energy density lt m in is provided. We do 
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not compute physical parameters for the sources unresolved 
in VLA images because no tight constraints on their val- 
ues can be derived. To compute the physical parameters 
of MRK 1239 and NGC 3783 we make use of the values 
at 1.6 GHz reported in Table 3. In the case of NGC 5506 
and NGC 7469 we compute th e physical parame t ers m aking 
use of the values rep orted in IMiddelberg et all l|2004 ) and 
lLonsdale et al.l (|2003l ), respectively. The brightness temper- 
ature is computed by: 

T B = W (cV 

where S(u) is the flux density at the frequency i/, # ma j and 
#mm are the source major and minor axis, k is the Boltzmann 
constant, and c the speed of light. For the m agnetic field 
and energy density we follow standard formulae l|Pacholczvkl 
1970), assuming equipartition condition between radiating 
particles and magnetic field. The equipartition magnetic 
field is 



H cq — y — 7m m in 

and the minimum energy density u m in is 

/ T \ 4/7 
u min = 1.04 ■ 10~ 23 ■ (-) 

where L is the synchrotron luminosity in Watt, and V the 
volume of the emitting region in pc 3 . We assume that the 
volume of the emitting region is a prolate spheroid 

V = "T^maj^min 

6 

where d ma j and d m i n are the source linear sizes. We assume a 
filling factor of unity, i.e. the source volume is fully and ho- 
mogeneously filled by relativistic plasma. Furthermore, pro- 
ton and electron energy densities are assumed to be equal. 
For the Seyfert nuclei with VLBI information, we obtain 
U min of the order of 10" 6 - 10~ 7 erg/cm 3 , H cq ~1.5 - 10 
mG and T B of about 10 7 - 10 8 K. These values are compa- 
rable _to_Uiosj^d«nwd_ta Seyfer t galaxies 
(e.g. lUlvestad fc Wilsonlll984l ; iKukula et afll 19991 ). 

5.2 The off- nucleus radio emission 

Of the seven Seyfert nuclei studied, four (NGC 1097, 
NGC 5506, NGC 7469 and NGC 7582) are surrounded by 
diffuse emission extend ing on kpc scale. In NGC 5506, 
IWehrle fc Morris! (119871 ) suggested that its extended 0.35- 
kpc halo may be due to either a radio plasma bubble 
expanding from the AGN, or a magnetically dominated 
coronal arc. Another possibility is that it is produced 
by free-free emission from the ionised gas heated by the 
AGN. However, the steep spectral index a = 0.9 makes the 
interpretation unlikely, as thermal emission would produce 
a flatter spectral index. Furthermore, the extended emission 
does not seem to be related to starburst activity, as no 
evidence for this is se eing in high spat i al resolution IR 
imag es up to 20 /im (|Prieto et al. 2009; R eunanen et al.l 
12009 ). 

In NGC 1097, NGC 7469 and NGC 7582, the extended 
emission is resolved in knots (Table 3), each having a 
counterpart star-forming region in high spatial resolution 



IR images l|Reunanen et al.ll2009l ). The analysis of the SED 
of the individual star-forming clumps in these galaxies, 
from radio to UV, is showing that these regions are young 
super-massive stellar clusters, an alogues to those seen in the 
starb urst nucleus of NGC 253 dFernandez-Ontiveros et all 
2009). In radio, the spectral index is rela tively steep, ~ 0.4 
- 0.9 in NGC 1097 (jHummel et al.lll987T l and in NGC 7469 
l|Alberdi et all 120061 ). and ~0.9 in region Rl of NGC 7582 
(no sufficient information is available for the other knots 
of NGC 7582 to derive the spectral index), suggesting a 
dominant non-thermal origin, most probably from super- 
novae events in these clusters. Not all the IR star-forming 
regions are detected in radio, suggesting different evolution 
phase of the clusters within the same star-forming ring 
(Fernandez-Ontiveros et al. 2009, in preparation). 



5.3 The missing flux density 

5.3.1 Observational limitations 

An interesting characteristic shown by the Seyfert nuclei 
with VLBI information, namely MRK 1239, NGC 3783, 
NGC 5506 and NGC 7469, is that their pc-scale flux den- 
sity is significantly lower than that measured on VLA 
images even in the presence of unresolved structures. In 
MRK 1239 and NGC 3783, the flux density derived from 
VLBA data accounts only for 20% of that measured on 
VLA images. This suggests that a significant fraction 
of the radio emission is not concentrated in a compact 
component at the centre, but it spreads over a region 
that is slightly larger than 0".l, i.e. the largest angular 
scale detectable by the VLBA at 1.6 GHz, and smaller 
than the resolution of VLA data at the same frequency, 
that is < ".3 for MRK 1239 fsee T abled) and < 0".5 for 
NGC 3783 (|Ulvestad fc Wilsorj[l983 ). A strong support to 
this interpretation comes from observations with different 
resolution of the nucleus of NGC 7469 . In MERLIN observa- 
tions at 1.6 GHz (|Alberdi et al.ll2006l ) the nucleus is slightly 
resolved on the east-west direction, giving an upper limit 
to its angular size of about 0".15. The east- west elongation 
is confirmed by pc- scale resolution VLBI observations 
l|Lonsdale et all 120031 ), where the nucleus is resolved in 
5 components spread over an area of about 0".17, i.e. 
comparable to that derived from MERLIN data. However, 
the flux density recovered from VLBI data accounts only 
for ~50% of that measured on M ERLIN images. A simila r 
result was found in NGC 5506 by IMiddelberg et all l|2004l) . 
where EVN observations of the nucleus could recover only 
43% of the flux density measured by MERLIN. In this 
case, an extended, low-surface brightness region, detected 
by MERLIN, is not seen by the EVN, as it is insensitive to 
structures larger than 35 mas at 18 cm, and 11 mas at 6 
cm. 

From these results we can argue that in the case of 
NGC 5506 and NGC 7469, the missing flux density on VLBI 
image is likely due to diffuse, steep-spectrum low-surface 
brightness emission that is undetectable by VLBI obser- 
vations either due to sensitivity limitation or because the 
structure is larger than the maximum scale detectable by 
the interferometer, or a combination of both. A similar 
result was found in the case of the Seyfert NGC 1068 where 
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IGallimore et al. I (120041 ) noted that deep VLBA observa- 
tions could detect a flux density almost 50% higher than 
that o f lower sensitivity VLBA observations by IRov et al. I 
(1998), but still the total flux density measured on images 
with lower spatial resolution could not be recovered. 



5.3.2 Steep- spectrum versus flat-spectrum Seyfert nuclei 

The fact that in Seyfert nuclei the pc-scale radio emission 
does not account for all the flux density measured at 
lower resolution, even in the pre sence of an unreso lved 
component, was already noted by Sadler et al ] |[l99l) by 
studying a sample of 22 nearby Seyferts with the PTI. 
They suggested that in Seyfert galaxies the radio emission 
is more diffuse and less centrally concentrated than in 
elli ptical rad i o gala xies. A similar behaviour was also found 
by lLal et aD ([2004) by comparing simultaneous VLA and 
VLBI observations of a sample of Seyfert galaxies. They 
found that in ~ 60% of the sources more than half of the 
VLA flux density of the unresolved component is missing 
in VLBI images. 

However, evidence of undetected flux density on pc-scale is 
not a characteristic common to all Seyfert nuclei. For exam- 
ple, it has been found that in flat-spectrum Seyfert nuclei 
essentially all the emission on arcse c ond scale is present on 
VLB I scales (jAnderson et al.l 120041 : lAnderson fc Ulvestad I 
2005), indicating that all the radio emission is concentrated 
in the central component. On the other hand, in Seyfert 
nuclei characterised by a steep spectrum, the flux density 
recovered on mas-scales is usually significantly lower than 
that expected from arcsecond-scale images, suggesting that 
the radio emission is not centrally concentrated, but it 
is diffuse on a larger region. A noticeable example of a 
steep-spectrum Seyfert nucleus where the flux density al- 
most disappears moving from arcsec to milli-arcsec scales is 
represented by NGC 4151, where only 8% of the flux density 
from the cent ral component could be recovered by VLB I 
observations (jPedlar et al. I 1 19931 : I Ulvestad et al] 1998), 
while a large fraction of the jet emission is undetected. This 
suggests that in Seyfert nuclei unresolved on arcsec scale 
but with different spectral index properties, the dominant 
radio emission may originate from different features: from 
the central core in flat-spectrum objects, from extended 
features, like jets, in steep-spectrum nuclei. 



5.3.3 Thermal or non-thermal origin? 

So far, the nature of the missing flux density has not been 
investigated in detail. This diffuse emission may be due 
to either thermal emission from ionised gas heated by the 
AGN, or by nuclear star forming regions, or non-thermal 
radiation from the AGN itself. 

To investigate a possible free-free origin of the missing flux 
we compute the electron density n e that the ionised gas 
must have in order to emit the "missing" flux density: 

n\ = 1.84 x 10 41 V2 DlS{u)V- x g^ (1) 



where Dl is the luminosity distance, T is the gas tempera- 
ture in units of 10 4 K, and gs the Gaunt factor. In Eq. [T]we 
assumed that electrons and protons have the s ame density. 
At ra dio frequency, g s ~ 17.7 + ln(T 3/2 /u) (e.g. IPacholczvkl 
Il97d ). We assume the gas temperature T ~ 10 4 — 10 6 K. 
The missing flux is calculated as the difference between the 
VLA and VLBI flux densities. In the case of MRK 1239 
and NGC 3783, the missing flux at 1.6 GHz is ~ 50 and 
18 mjy, respectively. If we consider these parameters in 
Eq. [TJ we obtain an upper limit to the electron density 
10 3 - 10 4 cm 3 , over a spherical volume of 0.1 arc- 
sec (i.e. 57 and 20 pc for MRK 1239 and NGC 3783 respec- 
tively) in diameter. This value is the largest angular size 
detectable from the VLBA, (i.e. the VLBA at 1.6 GHz is 
insensitive to structures larger than 0".l) and it represents, 
therefore, a lower limit to the source size. Such a dense gas 
would completely absorb the synchrotron emission arising 
from the embedded AGN: 

S ohs (v) = S (v)e-^ )l (2) 

where S b s and Sq(v) are the observed and the intrinsic flux 
density at the frequency v respectively, / is the width of the 
ionised re gion, and n(v) the absorptio n coefficient at the fre- 
quency v (|Rvbicki fc Lightman 1 19791 ). At radio wavelengths 

H(v) ~ 0.018T" 3/2 n^" 2 5ff . 

We consider that the pc-scale structure found in the VLBA 
images (Figs. 3t and[5| arises from the central AGN and the 
observed flux density is 10 mjy and 5 mjy for MRK 1239 
and NGC 3783 respectively. We assume that the AGN is 
located in the centre of the ionised region, i.e. I is equal to 
the radius of the ionised sphere (Fig. [9} that accounts for 
28.5 pc and 10 pc in the case of MRK 1239 and NGC 3783, 
respectively. With this parameters in Eq. we obtain an 
intrinsic flux density arising from the central AGN that 
is much higher than a few thousands of Jy. This suggests 
that the "missing" flux is not free-free from gas ionised by 
the central AGN, and a non-thermal AGN-related origin is 
more plausible. 

Another possibility may be that the ionised gas is related 
to an Hn region located in projection behind the AGN, or 
that the gas around the AGN is not ionised uniformly. In 
this way, our line of sight to the central AGN does not pass 
through the ionised gas, thus avoiding any absorption of 
the AGN radiation. 

However, a strong support to the non-thermal synchrotron 
origin of the diffuse emission comes from the comparison, 
when possible, between the spectral index distribution in 
low and spatial resolution images. For example, in the 
case of NGC 5506, the unresolved nucleus in VLA data 
has a steep spectru m with a mean valu e a ~ 0.8, while at 
pc-scale resolution ( Sa dler et al.l [l995;) the spectrum has 
a convex shape with the peak occurring around 3.5 GHz. 
This suggests that a significant fraction of steep-spectrum, 
low-surface brightness emission is present in the nuclear 
region, but it cannot be detected in VLBI observations 

likely due to observational l imitations. 

As in the case of NGC 4151 j|lllvestad et al.lll998l ). the steep 
emission may arise from a jet, that may be distorted and/or 
disrupted by the dense ambient medium. For example on 
arcsecond scale, NGC 4151 displays a compact core and tw o 
coUimated jets (|Pedlar et al. 1 1 19931 : IKukula et ai~l Il995l ). 
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Figure 9. The assumed geometry for the ionised gas heated by 
the AGN. The AGN is located in the centre of the ionised sphere 
with an angular diameter of 0.1 arcsec, that corresponds to the 
largest region detectable by the VLBA at 1.6 GHz. 



When observed with higher resolution, the compact core 
is resolved in several knots fo rming a .jet which i s not 
aligned with the arcsec-scale jet l|Ulvestad et al.l ll998). The 
misalignment often found between the radi o jets on milli- 
arcsecond scale and those on larger scale l|Ulvestad et al.l 
1 19981 ) may be an evidence of the interaction between 
the radio jet and the environment. In this scenario it is 
also possible that, for some reason, flat-spectrum Seyfert 
nuclei are not able to develop a radio jet, and this causing 
that all their emission is produced in the central compact 
component. 



6 CONCLUSIONS 

We presented the analysis of multi-frequency archival VLA 
and VLBA data of 7 close and bright Seyfert nuclei. The 
conclusions we can draw from this investigation are: 



• At VLA resolution, FWHM ~ 0".l, the nucleus of 
the Seyfert galaxies is unresolved, with the only exception 
of MGC-5-23-16 and NGC7469 which show a core-jet 
structure. At VLBA resolution, equivalent to a few parsecs, 
the nucleus of MRK 1239 is resolved in two components, 
while in NGC 3783 it is still unresolved. 

• The Seyfert galaxies in this study with known circum- 
nuclear star-forming regio ns in the IR, namely NGC 1097, 
NGC 7469 and NGC 7582 (|Reunanen et al.ll2009l ) present a 
radio counterpart for a few (^10%) of these regions. Most 
are not detected in radio, indicating that those detected 
may be characterised by higher supernovae rate. Their steep 
radio spectral index is in line with this idea. Conversely, 
none of the other Seyfert nuclei present any circumnuclear 



radio emission. The only exception is NGC 5506 that shows 
a radio halo surrounding the nucleus. 

• A comparison between arcsecond and milli-arsecond 
resolution in MRK 1239 and NGC 3783, pointed out that 
almost 80% of the radio flux density detected in VLA 
observations is not recovered at pc-scale resolution. This 
suggests the presence of a diffuse component on scales of 
a few tens of parsecs, undetected with the VLBA. Similar 
situation is found in NGC 1068, NGC 4151, NGC 5506 and 
NGC 7469. The nature of this "missing" flux components 
is likely due to synchrotron AGN-related emission. This 
difference between the flux density measured on arcsecond 
and milli- arcsecond resolution images is not found in the 
case of elliptical radio galaxies, but it appears to be a 
common phenomenon in Seyfert galaxies with a steep 
spectrum, mostly hosted in spirals. If of synchrotron origin, 
this emission may be spilt off from a jet that may get 
easily distorted and/or disrupted by the dense interstellar 
medium in the nucleus of spirals. 

• A comparison between Seyfert nuclei with different 
spectral properties points out that in flat-spectrum nuclei, 
almost all the flux density is recovered on milli- arcsecond 
scale. This indicates that in flat-spectrum objects the ra- 
dio emission is essentially concentrated in the compact core, 
without evidence of jet-like structure even on milli-arcsecond 
scale, while in steep-spectrum objects a significant fraction 
of the radio emission arises from low-surface brightness, ex- 
tended features. 
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